Exploiting lattice potentials for sorting chiral particles 
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Several ways are demonstrated of how periodic potentials can be exploited for sorting molecules 
or other small objects which only differ by their chirality. With the help of a static bias force, the 
two chiral partners can be made to move along orthogonal directions. Time-periodic external forces 
even lead to motion into exactly opposite directions. 

PACS numbers: 05.40.-a 05.60.-k,05.45.-a 



Chiral particles are extended objects which are non- 
supcrposable with their mirror image. So-called enan- 
tiomcrs, i.e. chemically identical molecular species with 
opposite chirality play a crucial role in Chemistry, Bi- 
ology, and Medicine due to the omnipresence of chiral 
molecules in living organisms but with only one of the 
two chiral partners actually being present. Accordingly, 
cnantiomers in drugs, pesticides etc. have very differ- 
ent effects on an organism and thus their separation is 
of great importance. Established methods of separating 
cnantiomers mostly exploit some kind of chiral selector 
[!||, i.e. some materials, structures, or ancillary molecules 
which themselves exhibit an intrinsic chirality. Their 
main disadvantage is that essentially every enantiomer 
species requires a different selector. Therefore, several 
alternative concepts have recently been put forward. A 
first promising direction proposes to utilize appropriate 
microfluidic flows, such as vortices @ or shear flows 
Second, photoinduccd separation by means of suitably 
chosen electromagnetic fields has been theoretically pre- 
dicted in Ref . Q . A third approach to exploit a structure 
without an intrinsic chirality is due to de Gennes , pre- 
dicting qualitatively that, according to Curie's principle 
0, small chiral crystalls should slide down an inclined 
plane along directions which slightly differ for the two chi- 
ral partners, provided thermal noise is negligible. Here, 
we further pursue this approach, showing that with the 
help of periodic potentials the two chiral partners even 
can be made to move into opposite directions, with re- 
markable persistence against thermal noise. 

Apart from "true" (bio-) molecular cnantiomers, we 
also have in mind chiral nano- and micro-particles, e.g. 
helically shaped nonmotile bacteria Q and artificial flag- 
ellae §, carbon nanotubes, chiral colloidal clusters @, or 
ferromagnetic nano-propellers The periodic poten- 
tials we are proposing to utilize for sorting those chiral 
particles may be realized e.g . by means of crystal sur- 
faces [ll|, optical lattices [l2f . periodic micro- and nano- 
structures [IH, or magnetic bubble lattices flij ]. 



FIG. 1: (a) Solid: Triangle, specified by a, 6, c. Dashed: Its 
chiral partner, (b) and (c): Direction •& and modulus v of the 
net velocity v — e$v versus the direction a of a static bias A in 
([5]) by numerically solving (J2j|- ((6j) with parameters as specified 
below ©. Solid: kT = 0.08, A = 3.6, a = 0.4, b = 0.23, 
c = 0.34 (proportional to the solid triangle in (a)). Dashed: 
Same but for the chiral partner (dashed in (a)). Dashed- 
dotted and dotted: Same but for A — 10. 



planar, chiral "molecule", consisting of three identical, 
rigidly coupled "atoms" or other small objects with bro- 
ken mirror symmetry, see Fig. la. All basic effects and 
mechanisms are recovered in three dimensions and also 
for more general lattices, but are much more cumber- 
some to visualize and explain. This very general validity 
of our main results will be exemplified for various other 
chiral "molecules" in the end. We thus consider the two- 
dimensional dynamics 



rnAit) = -jiXi{t) + F(xi(t), t) + /i + &(t) 



(1) 



Most of us are not very accustomed to think in terms 
of chiral symmetry and symmetry-breaking, especially in 
combination with the crystal symmetries of a periodic 
potential. For this reason only, we mainly focus on the 
simplest possible setup 0, namely the two-dimensional 
dynamics in a square lattice potential of a "minimal" 



Dots indicate time derivatives, xi = e.\x^\ + e^a^a arc 
the "atom positions" (i = 1,2,..., A) in Cartesian coor- 
dinates e„ (u = 1,2), mi their mass, and ji their dissi- 
pation coefficient, e.g. due to an ambient fluid. In par- 
ticular, for the triangular particles (Fig. la) we have 
N = 3 and z-independent 7, and rrii. The force field 
F(x,t) is partly due to a "lattice potential" (see be- 
low) and partly due to an externally applied driving, 
typically via electrophoresis. Under these conditions, 
hydrodynamic interactions are screened [Tij and there- 
fore safely negligible 0, The internal constraining 
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forces, maintaining e.g. the triangular shape in Fig. 
la, are represented by /j, and thermal fluctuations are 
modelled as usual by unbiased Gaussian white noise 
£i(t) = ei£i,i(t) + (32^,2 (t), satisfying the fluctuation dis- 
sipation relation (£,i,fj,(s)t;j t v(t)) = 2jikTSijS fJll/ S(s — t) 
with T the ambient temperature and k Boltzmann's con- 
stant. The position of the rigid "molecule" is conve- 
niently specified by the so-called center of friction Q 

X := ^ i= i SfcLi 7*j ana - ^ s orientation by the an- 
gle cj> between the e*i axis and the position of "atom 1" 
relative to the center of friction: (f> := <(ei,afi — X). 
Rewriting (JTJ in terms of X and to get rid of the con- 
straining forces fi is a basic mechanics exercise. Further, 
for the very small objects we have in mind, inertia effects 

are nee hgibie mEa, y ieldin g 



X(i) 
0(t) 



2ji=i 7i 2/j 

, m = O(<t>(t))yi(0) . (4) 



(2) 



+ (3) 



In ([3J, vectors are temporally embedded into R 3 with 
standard scalar and vector products ■ and x . In (QJ, O(0) 
is a rotation matrix with elements On — O22 = cos <f> and 
O12 = — O21 = — sincf). Thus, yi(i) are the particle pos- 
tions relative to the center of friction with convention 
</>(0) = and with t- independent modulus %ji := \yi(t)\. 
Finally, and C/>M are independent Gaussian white 
noises with (C M (s)G(*)} = 2kT8 llv 5{s - t)/Y,li and 

(C*(*)C*(*)> = 2fcT*( fl -t)/E7<W?- 
As already said, the force field consists of two parts, 



F(x,t) = e a A(t) - VU(x) , 



(5) 



namely a spatially homogeneous, externally applied force 
along the direction e a := e± cos a + e*2 sin a and a Gaus- 
sian square lattice potential with period L: 
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For this potential with u > and u < as well as for 
various other potentials we always found similar results. 
Focusing on u > from now on, the natural energy scale 
is the potential barrier A(7 := U(eiL/2) — U(0) separat- 
ing adjacent potential wells. We henceforth adopt time, 
energy, and length units so that mini7i = 1, AU = 1, 
L = 1, and focus on a = L/4 (Til ]. 

The quantity of central interest is the net velocity v = 

e#v, obtained by averaging X{t) over time. Obviously, 
rotating the force field ([5]) by 90° leaves the potential © 
invariant and entails a rotation of v by 90°. Hence, it is 
sufficient to focus on a G [0°,90°]. Likewise, one readily 
sees that A(t) n- —A(t) implies v n> —v. 

We first consider ^independent A, i.e. the force field 
([5]) derives from a tilted periodic potential. For A = 



symmetry implies v = 0. For 4^0 the salient point is to 
realize that there exists no symmetry argument why two 
"molecules" of opposite chirality should travel down the 
tilted periodic potential with identical velocities v. Fol- 
lowing de Gennes @ , we thus can invoke Curie's principle 
to conclude @, [H| that generically (i.e. up to parameter 
sets of measure zero) the velocities will indeed be dif- 
ferent. In other words, (practically) any tilted periodic 
potential can separate chiral partners via their velocities. 
The main remaining problem pinpointed by De Gennes 
is the quantitative efficiency of the effect. 

Fig. 1 provides those quantitative details in a typical 
case. We see that the velocities v of the two chiral part- 
ners are indeed disappointingly similar, except around 
a = 45°. The explanation is as follows: For small ther- 
mal energies kT and small bias A, the particles travel 
extremely slowly by thermally activated hopping from 
one local minimum of the tilted periodic potential to the 
next. For any given orientation a there exists a critical 
tilt A in ([5]) at which certain local minima disappear by 
annihilation (collision) with saddle points, giving rise to 
"running solutions". For kT = (deterministic limit), 
these solutions travel either parallel to e\ or to e*2, and 
for small kT > still almost so. Roughly speaking, the 
direction "closer" to that of the static bias e a A is pre- 
ferred, but due to the broken mirror symmetry, the di- 
rection actually switches already at some a < 45° for 
one chiral partner and symmetrically at a > 45° for the 
other (solid and dashed in Fig. 1). Since these consid- 
erations do not depend on any details of the model we 
can conclude that a separation by (almost) 90° is generic 
for a = 45°, small kT, and A close to criticality. Upon 
further increasing A, the deterministic running solutions 
speed up and bifurcate into new ones, "locked" [l2|, [I?} 
along directions of the form ne\ + me,2 with increasingly 
large integers n and m (dashed-dotted and dotted in Fig. 
1) and with v — > e a AN/ ^7i f° r A — > 00. Likewise, for 
finite kT the deterministically "sharp" bifurcations get 
washed out (Fig. 1) and v — > e a AN/ ^2^fi for kT — > 00. 

Thus focusing on a = 45°, the dependence of the ve- 
locity v on the bias A is shown in Fig. 2 (a,b). For sym- 
metry reasons, the velocities of the two chiral partners 
are now equal in modulus and symmetric about a = 45° 
(see also Fig. 1). Remarkably enough, for some A- values, 
one triangle moves (practically) parallel to e\ (Fig. 2c) 
and thus its chiral partner parallel to e.2 (not shown in 
Fig. 2), while for some different A- values it is exactly the 
other way round (Fig. 2d). In other words, one and the 
same triangle may move along orthogonal directions for 
two different A-values. 

Turning to periodic A(t) in ([5]), our so far findings quite 
naturally suggest the following idea: We select a = 45° 
and two static bias values A\ and A2 with velocities 
v\ = v\<=\ and V2 = (e.g. A\ = 4 and A2 — 6 

for the solid lines (kT = 0.014) in Fig. 2a, b) and exploit 
that the signs of v\ and can be arbitrarily chosen by 
adjusting the signs of A\ and A2 (recall that A —A 
implies v i— > —v). If we now construct a time-periodic 
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FIG. 2: (a) and (b): Direction i? and modulus t> of the net ve- 
locity v = e$v for the solid triangle from Fig. la versus static 
bias A with fixed direction a = 45° . Shown are numerical so- 
lutions of ©-© for kT = 0.014 (solid), kT = 0.16 (dashed), 
and kT = 0.32 (dotted). Other parameters as in Fig. 1. (c) 
and (d): Illustration of the triangle's motion for A = 4 (c) 
and A — 6 (d) at kT — 0.014. The static bias is indicated 
by A := e a A and the periodic potential as "shaded back- 
ground". The triangle motion is shown for a time-span of 
about 14.5 time-units in (c) and about 8.5 in (d), and then 
continues periodically up to noise effects (not shown). 



A(t) which takes the value A\ during a fraction p G [0, 1] 
of its total period r and the value A 2 during the rest of 
the period, the resulting time averaged velocity will be 
v = pvi + (1 — p)v2, provided r is so large that transient 
effects after each jump of A(t) are negligible. Therefore, 
the molecule can be steered into any direction on the two- 
dimensional plane by varying p and adapting the signs of 
A\ t 2- In particular, we will encounter a situation where 
v is orthogonal to the force direction e a . E.g. from the 
solid lines in Fig. 2 we can read off that the triangle will 
move with such a velocity v _L e a if we choose A\ = —4, 
A2 = 6, and p ps 2/3 to account for the difference in 
modulus of the corresponding velocities v\ ps — 1.2ei 
and V2 Ps 2.4 e%. The net velocity v of the chiral part- 
ner follows from the above mentioned symmetry about 
a = 45°: This symmetry applies to both v\ and v*2 sepa- 
rately, and hence also to v = pv\ + (1 — p)v2- Altogether, 
the two chiral partners can thus be forced to move into 
exactly opposite directions. Deviations due to the so far 
neglected transient effects after each jump of A(t) are - 
at least for not too small r-values - small and thus can 
be compensated by adjusting p and/or A\^- 

Fig. 3 shows that these ideas indeed work out in prac- 



0.5 




0.7 



0.65 I 





1 1 


1 

00 - 




. (a) r = 


12 




^\ T — 


6 






3 




- - - r = 


1 


1 

— 


1 h— 

1 1 


1 



0.02 0.04 0.08 



kT 



0.32 0.64 





62 _ 




e*2 _ 




e*2 ^ 


(c) 

1 — 


\ 


(d) 

h 


X 

H 


(e) 
1 


— \\ 



FIG. 3: (a): Absolute velocity v versus thermal energy kT 
(logarithmic scale) for the same system as in Fig. 1 but with 
a r-periodic driving e a A(t) with a = 45° and A(t) taking 
the value A\ = —4 during a fraction p of the period r and 
the value A2 = 6 during the rest of the period, (b): The 
corresponding p- values, adjusted as described in the main text 
so that the two chiral partners move into exactly opposite 
directions, (c) Typical single-particle trajectories X (t) for the 
two chiral partners (light red and dark blue) with t G [0, 100], 
X(0) = 0, kT = 0.02, p = 0.66, and r = 6. Other parameters 
as in (a). The bar indicates 50 lattice periods and the double 
arrow the periodic driving (d) Same but for a much 

larger thermal energy kT = 0.32 and t G [0, 6000], p = 0.69. 
(e) Same as in (c) but for a very different triangel with a — 2.5, 
b = 2.2, c = 1.1 (cf. Fig. 1), and a driving with Ai = -7, 
A 2 = U,p = 0.8 pi]. 



tice, and in fact down to surprisingly small time-periods 
t and up to remarkably large thermal energies kT. Note 
that while the velocities in Fig. 3a arc long-time aver- 
ages, Figs. 3c-e exemplify single-particle trajectories of 
moderate duration. Hence the thermal noise still leads 
to quite notable random fluctuations of each trajectory 
X(i) around the average behavior, especially in Fig. 3d. 
Only in Figs. 3c, e we still can see the expected "steps" 
of X(t) at jumps of A(t). 

Our above recipe for tailoring transport directions can 
be readily extended to arbitrary velocities v\ and V2, pro- 
vided they are not parallel to each other: Then, as be- 
fore, v = pv\ + (1 — p)v2 can be made to point along any 
direction by properly choosing p and the signs of A\^- 
Intuitively and in view of Fig. 2, it is quite clear that 
generically one will always be able to find two bias values 
A\ and A 2 with non-parallel velocities v\ and V2 ■ We thus 
can conclude that chiral partners can (practically) always 
be made to move into opposite directions by means of a 
suitably tailored periodic driving force. Fig. 3e exempli- 
fies this generalized theoretical scheme for comparatively 
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FIG. 4: (a)-(c): Further examples of chiral "molecules". 
Dots indicate the constituting "atoms", lines their rigid cou- 
pling, and the "shaded background" the periodic potential 
©. Adopting units as specified below (|6]), ji = 1 for all 
"atoms" except for the larger dots in b and c, represent- 
ing ji — 3. (d) Typical single-particle trajectories X(t) for 
the molecule from (a) and its chiral partner (light red and 
dark blue), obtained by numerically simulating ©-© with 
* € [0, 100], X(0) = 0, and kT = 0.02. Parameters of the peri- 
odic driving (see main text): a = 45°, At = 7.8, A 2 — —11.7, 
r = 6, p — 0.71. The bar indicates 50 lattice periods and 
the double arrow the periodic driving, (e): Same but for 
the "molecule" from (b) and Ai = 7, A 2 = -10.5, p = 0.78 
[Isl ]. (f): Same but for the "molecule" from (c) and A\ — 3, 
A 2 = -6, p = 0.85 H. 



"large" triangular particles and Fig. 4 for a representa- 
tive selection of more general chiral "molecules". Gen- 
eralizations involving more than two "static velocities" 
Vi and the concomitant optimization problems point into 
interesting directions for future research. 

In conclusion, periodic potentials can act as very effec- 
tive and versatile selectors for sorting small objects which 
only differ by their chirality. Static bias forces make the 
two chiral partners move into directions which differ by 
up to 90° (Figs. 1,2). Appropriately chosen time-periodic 
forces even lead to motion into exactly opposite direc- 
tions (Figs. 3,4). A major advantage compared to many 
other separation concepts [l[ is that one and the same pe- 
riodic potential may act as an efficient selector for quite 
different chiral particle species by suitably adapting the 
time-periodic driving force. Furthermore, the separation 
mechanisms are remarkably robust against thermal noise. 
The basic symmetry breaking conditions at the origin of 
all these effects arc gcnerically satisfied for much more 
general systems than in (U)-©, including three spatial 
dimensions, finite inertia effects, other chiral objects and 
crystal potentials. An experimental proof of principle 
for chiral micro-particles [8l4l0j moving in a periodically 
structured microfluidic device [13J is presently under con- 
struction in the Ansclmctti lab at Bielefeld University. 
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